Hyperfine interactions in silicon quantum dots 
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A fundamental interaction for electrons is their hyperfine interaction (HFI) with nuclear spins. 
HFI is well characterized in free atoms and molecules, and is crucial for purposes from chemical 
identification of atoms to trapped ion quantum computing. However, electron wave functions near 
atomic sites, therefore HFI, are often not accurately known in solids. Here we introduce an all- 
electron calculation for conduction electrons in silicon and obtain reliable information on HFI. We 
verify the outstanding quantum spin coherence in Si, which is critical for fault-tolerant solid state 
quantum computing. 
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I. INTRODUCTION 

The current interest in silicon in the context of spin- 
tronics and spin quantum computation^^ arise naturally 
from its perceived technological and fundamental ad- 
vantages. Technological advantages included scalability 
and unique material processing capabilities, developed 
primarily by the well established and, already for half- 
a-century dominant Si microelectronics industry. The 
fundamental advantage, which is the key motivation for 
our present work, is the very low abundance of finite- 
spin isotope ( 29 Si) in natural Si, leading to very long 
electronic spin coherence times and thus enabling fault- 
tolerant spin qubit operations. Recent ESR experiments 
on phosphorus donor electrons in isotopically purified 
Si show that electron spin coherence times are longer 
than 60 m4, and there have been tremendous recent ex- 
perimental progress in the study of the Si:P system.— ~— 
However, fabricating P donor arrays with the necessary 
atomic-scale precision^ is a difficult challenge still to be 
overcomo^— . An alternative to the Si:P system uses 
electron spins in gate-confined Si quantum dots (QD) as 
qubits. There have been extensive experimental stud- 
ies of few-electron spin dynamics in gated semiconductor 
(in particular GaAs) QDs for the past decade^, while 



recent experiment! 
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have shown that single-electron 



QDs are well within reach in Si as welli£. 

In Si:P, quantum coherence of donor confined electron 
spins is limited by their well-known hyperfine interactions 
(HFI) with the random 29 Si nuclei in the environment, 
whose dipolar interaction induced nuclear spin dynamics 
and the subsequent electron spin spectral diffusion is the 
main source of donor electron spin decoherence^i^— 
Clearly, accurate knowledge of HFI would also be re- 
quired to evaluate spin decoherence for an electron con- 
fined in a QD and to assess the magnitude of the Over- 
hauser field, i.e. the net field felt by the electron spin due 
to all the 29 Si nuclear spins in a Si QD. The strength of 



HFI in a QD with smooth confinement could be directly 
derived from the HFI of an extra conduction electron in 
bulk Si. In spite of its considerable importance, however, 
these HFI parameters have never before been evaluated 
for Si conduction electrons, and our work fills this infor- 
mation gap. 

Motivated by spintronic and quantum computing con- 
siderations, we develop an ab initio study of the HFI 
parameters between Si conduction electrons and 29 Si nu- 
clei. Our results are directly relevant to the prospective 
application of Si as the material of choice in these ap- 
plications. In addition, they arc in good agreement with 
the existing experimental data related to HFI^ 3 - - — Our 
results provide the basis for a local spin density approxi- 
mation (LSDA) for the HFI, similar to LDA in electronic 
calculations. The LSDA is ideal for treating electrons in 
relatively extended QD states in Si, as illustrated below. 

We emphasize that the calculation of the HFI for a 
single extra electron in an unperturbed system, like bulk 
Si, introduces fundamental methodological difficulties as 
compared to the spatially localized electron bound to 
shallow impurities in semiconductor s 1 ^ 20 ' 26 " — due to the 
extended nature of the bulk electron. In general donor 
electron wave functions cannot be expressed simply as an 
envelope function multiplied by the Bloch states of the 
conduction electrons^ Therefore, even though calcula- 
tion of HFI at an impurity site or at neighbouring sites for 
shallow impurities in semiconductors is straightforward 
nowadays ) 20 ' 26 ~ — we cannot easily extract information 
about the underlying Bloch functions, which determine 
the conduction electron HFI. 

The paper is organized as follows. In Section II we 
introduce the theoretical method with which we perform 
the hyperfine calculation for bulk Si. In Section III we 
show the results on conduction electron hyperfine inter- 
action in Si and compare with existing experimental mea- 
surements. In Section IV we calculate the hyperfine in- 
teraction for confined electron in a Si quantum dot, and 



compare our results with that in GaAs quantum dots. 
Lastly, in Section V we draw our conclusions. 



II. THEORY AND FORMALISM 



The HFI hamiltonian has the general form 



#HF = J • A ■ S , 



(1) 



where S and I are the electron and nuclear spin op- 
erators, and A is the HFI tensor. The tensor compo- 
nents for a nucleus at Rj can be written in terms of an 
isotropic term a, known as Fermi contact term, and of an 
anisotropic dipolar traceless tensor bij : Aij = aSij + 6, y . 
The labels i and j refer to coordinates x, y, and z taken 
here along the conventional fee cubic cell edges. We con- 
sider Si with one extra electron at a fixed conduction 
band minimum k{. The system has axial symmetry with 
respect to the i-axis, so that the traceless anisotropic ten- 
sor is diagonal and expressed as b^ = diag{2b, —6, — b}. 
The hyperfme tensor A is thus completely defined by 
the scalar hyperfme parameters a and 6, as described in 
Ref. |30j. Explicit expressions for these parameters^ are 
determined by the electron spin density, 



Ps(r) =P\{r)- Pi{r) 



(2) 



In particular the Fermi contact interaction is propor- 
tional to the spin density at the nuclear site, ps(Ri)- 
From a and b, given in energy units, one can obtain the 
equivalent magnetic field created by the nuclear spin act- 
ing on an electron spin. Reliable calculation of the HFI 
parameters require precise values of ps{f) for all r, in 
particular in the vicinity r ss Rj. 

Our calculations, performed within the Density Func- 
tional Theory (DFT) framework,— involve the full- 
potential linearized augmented plane wave method (FP- 
LAPW)^ as embodied in the WIEN2k packaged This 
state-of-the-art all-electron methodology includes spin 
polarization of the core and valence states, spin-orbit 
coupling, and relativistic effects. No shape assumption 
(e.g. the usual spherical constrain) to the potential is 
involved. The generalized gradient approximation of 
Perdew, Burke and Ernzerhof (GGA-PBE) is used for 
the calculation of the exchange-correlation potential^ 
Calculations of HFI of negatively charged bulk Si are 
performed for several Si supercell sizes, ranging from N 
= 8 to 64 atoms, each with a single extra electron. In 
practice this corresponds to different average electronic 
densities p^ — 2e/(NQpc), where ttpc is the volume 
of the Si primitive cell (PC), which contains two atoms. 
A uniform positive jellium background with density \pn\ 
is included to cancel the long range multipole interac- 
tions of charged supercells^ Grids from 8 to 125 k-points 
are used to sample the irreducible wedge of the Brillouin 
zone^ according to the supercell size. The calculation is 
performed self-consistently until convergence on the to- 
tal energy and total charge is achieved. This approach 



realistically describes the electronic properties, and has 
previously been used to calculate HFI of impurities in 
semiconductors ! 27 ' 28 Our calculations give a Si lattice pa- 
rameter of 5.47 A [experimentally measured value (exp.) 
is 5.43 A], bulk modulus B = 90 GPa (exp. 98.8 GPa), 
and cohesion energy E = 4.55 eV (exp. 4.63 cV). 

We follow a two-step procedure to compute the hy- 
perfme parameters at a given Si site for a supercell with 
an extra electron: 1) The converged charge density is 
obtained from a standard spin-polarized self-consistent 
calculation. 2) This charge density is used to get the 
potential energy of the Kohn-Sham equation, and the 
Hamiltonian is written assuming a fixed ki correspond- 
ing to one of the six conduction band minima. In this way 
we constrain the extra electron to a fixed k-point, in anal- 
ogy with the case of impurities in semiconductors j2£t— 
From the sum over all occupied bands n of |0^(r)| 2 and 
of |<^ (r)| 2 for each spin orientation a and a (where (j)^'' 7 
are the Kohn-Sham eigenstates), we obtain the net s pin 
density [Eq.@], as described in the appendix of Ref. [27|. 
Note that for neutral bulk Si ps = 0, so the calculated 
Ps{r) is entirely due to the added electron. This proce- 
dure provides as good quantitative estimates of spin den- 
sity as possible within available theoretical treatments. 



III. HYPERFINE INTERACTIONS FOR A 
CONDUCTION ELECTRON IN BULK SILICON 

Figure [1] shows the spin density in the (100) crystal 
plane calculated with an 8-atom supercell for negatively 
charged Si, where the extra electron is kept at k = k z at 
the conduction band edge. Note the anisotropy of the dis- 
tribution, elongated in the axial z direction. Comparison 
of this spin density with the charge density in the same 
plane for the conduction band k z pseudo-wavefunction, 
given in Fig. 1(b) of Ref. [38|, illustrates the similarity of 
the two distributions, except for regions just around the 
atomic sites. This is because results in Ref. |_3S| were ob- 
tained within the pseudopotential approach, where the 
core region potentials arc very different from the actual 
ones. The same similarity is obtained in any other plane. 
These considerations, and the fact that /Cg(r) is nonzero 
with the added electron, allows us to identify ps(r) in 
our all-electron calculation with the charge distribution 
of the extra electron, which in a one-electron scheme cor- 
responds to I^Pfe^r)! 2 , the absolute value of the Bloch 
function squared. 

Figure[2]summarizes the calculated 29 Si HFI as a func- 
tion of 1/N in the negatively charged supercells. Clearly, 
the HFI parameters a and b scale linearly with the elec- 
tronic density. For natural Si, the isotropic parameter a 
is of the order of a couple of mT, while the anisotropic 
parameter satisfies b/a ~ 3% for all N. We have also 
calculated the HFI parameters a and b for an N = 8 
supercell using the more recently developed GGA-WC 
exchange-correlation functional^ The resulting HFI pa- 
rameters show only minor changes: a (b) increases (de- 
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FIG. 1: (Color online) Electron spin density in a (100) plane 
for bulk Si with an extra electron in the conduction band 
minimum at fc z . The vertical axis is z and the color scheme 
runs from red (high density) to purple (low density), following 
the rainbow sequence. The circular high density spots are the 
Si atomic sites. 
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FIG. 3: (Color online) Calculated r\ for different supercell sizes 
(red squares). The average value and standard deviation are 
given by the lowest (red) asterisk and respective error bar. 
Experimental results are given by the upper (blue) horizontal 
dashed lines, following the average values of Refs.|23| and X2A, 
labeled S&W (1956) and Wilson (1964), respectively. The 
error bar for Wilson was obtained as an upper limit of the 
error estimated from Ref . [241. where 77 si is given in an expres- 
sion related to 77 gc The horizontal position of the asterisks 
is arbitrary, as they represent average values. The estimated 
value of 77 > 300 in Ref. [25| is off the scale here. 



Si, which contains no unsaturated bonds, the systematic 
error should be smaller than 20%. 

Within a mono-electronic framework, the relative 
weight of the conduction band wavefunction at a Si 
atomic site Ri = may be quantified by a parameter (77) 
defined below i^ 3 - It is widely adopted in experiments and 
in effective mass calculations. The original definition in 
a one-electron context and the analogous expression for 
the all electron framework here are, respectively 



FIG. 2: (Color online) Squares give the calculated hyperfine 
parameters versus 1/N, where N is the number of Si atoms 
in the supercell. The asterisks give the interpolated values 
expected for natural Si (~ 5% of 29 Si), with a fa 1.9 mT 
and b fa 0.06 mT. Solid straight lines are ajv = (37 /N) mT 
and fejv = (1.1/JV) mT, and give good fits to the data points. 
These results are the basis of a local spin density approxima- 
tion used in our quantum dot calculations, similar to LDA in 
electronic calculations. 



creases) by less than 5% (8%) as compared to the values 
shown in Fig. [5] Indeed, previous studies comparing dif- 
ferent exchange-correlation functionals showed that DFT 
methods are generally quite consistent with each other, 
and have at most about 20% systematic error compared 
to experimental observations^ We expect that for bulk 
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l*(Q)l 2 



^ VAE 



Ps(0) 
(ps)av 



(3) 



where in the denominators (Q)av is the average of Q over 
the normalization volume. We omit the subscripts "IE" 
and "AE" below. Ref. [H inferred 77 = 186 ± 18 for 29 Si 
from NMR experiments. Later work2^ identified an error 
in Ref. [23| and gave a corrected value of 77 = 178 ± 31. 
Ref. [25| deduced 77 by extrapolating Overhauser shift data 
and obtained roughly twice the value reported in Ref. [23|, 
with 77 > 300. Our calculation constitutes the first the- 
oretical estimate for 77 and give 77 = 159.4 ± 4.5. The 
results are presented and compared to the experimental 
estimate s 23 ' 24 in Fig. [31 



IV. HYPERFINE INTERACTIONS FOR A 

CONDUCTION ELECTRON CONFINED IN A 

QUANTUM DOT 

We now consider an electron confined in a QD near a 
[001] interface with a barrier material. In this case the 
±fc 2 valleys are lower in energy and make up the elec- 
tronic ground state. By symmetry, ps is the same for k z 
and — k z , so the results here apply to general superposi- 
tion states involving ±k z . In the envelope function ap- 
proach, the electron wave function in the QD is given by a 
bulk state modulated by a slowly varying envelope func- 
tion F(r), with [ ,, \F\ 2 dr = 1. In the all-electron 

v ' ' Jail space I I 

scheme, the envelope modulation is |-F| 2 , since l^l 2 is as- 
sociated with ps, while the bulk contribution to the spin 
density is normalized within a PC, so that the overall 
spin density in a QD (normalized in all space) is 



p Q s D (r) = \F(r)\^ PC p P s c {r). 



(4) 



In essence, our ab initio calculation for a supercell of vol- 
ume V is equivalent to having a uniform envelope func- 
tion F = 1/VV and a bulk spin density normalized in a 
PC. In a QD, F is non-uniform and extends over tens of 
nanometers. However, this change does not affect the cal- 
culation for the contact part of the interaction, because 
both the envelope function and the total spin density are 
normalized over the QD. 

The contact HFI field due to a nuclear spin at Ri in a 
QD is 



a(Rj) = a 2 |F(Rj)| 2 fi PC . 



(5) 



where a 2 is extrapolated from our supercell calculations 
(see Fig. [5]) to a PC (a supercell with TV = 2). Since we 
fit our calculations as ajv ~ (37 /N) mT, a 2 = 18.5 mT. 
The a(Rj) thus calculated is the interaction strength of 
the electron with a particular nucleus at R/ . If we have 
a QD where every nucleus has a spin, with the nuclear 
spins arranged in a fully polarized state, the total HFI A 
is a sum over the whole QD: 



A = J2 a ( Rl ) = 2a2 

R^QD 



(0) 



where the factor 2 accounts for contributions from the 
two nuclear sites per PC, and, to an excellent approxi- 
mation, the envelope is taken to be constant within each 
elementary PC. The total HFI A corresponds to the mag- 
netic field acting on the electron if all nuclear spins in the 
QD are polarized. 

In Table Q] we show estimates of the relevant Over- 
hauser fields 7 and/or energy scales in Si and GaAs QDs. 
The random Overhauser field 8 A for a high temperature 
nuclear reservoir is 5 A = A/y/Ms, with Ms giving the 
total number of finite spin nuclei. The T 2 time is the 
electron spin dephasing time, given by T 2 * = h/SA. The 
first row of the table shows typical data for a GaAs^ 9 - 
QD that contains about 10 6 nuclei, and is given as a 



host 


Mr 


Ms 


A 


SA 


T 2 * 


GaAs 


10 6 


10 6 


92 peV (3.6 T) 


92 neV 


7.2 ns 


Natural Si 


10 5 


5000 


210 neV (1.85 mT) 


3.0 neV 


0.22 ps 


100% 29 Si 


10 5 


10 5 


4.3 peV (37 mT) 


13.6 neV 


49 ns 


0.01% 29 Si 


10 5 


10 


0.43 neV (3.7 pT) 


0.136 neV 


4.9 ps 



TABLE I: Order of magnitude estimates of HFI that an elec- 
tron can experience in a Si QD. The columns give the total 
number Mt of nuclei in the dot, the number Ms of non-zero 
spin, the maximum A and random 8A — A/\/Ms Overhauser 
fields, when the nuclear spins are all aligned or in a random 
configuration, and the corresponding T 2 * dephasing time. The 
estimates for GaAs are given for comparison. 



benchmark for comparison. The 2nd row gives the data 
for a typical natural Si QD with 10 5 nuclei (generally Si 
QDs are more strongly confined than GaAs QDs). With 
the three orders of magnitude difference in A in GaAs 
and natural Si, the random fields and T 2 * arc different 
by about two orders of magnitude. For an isotopically 
enriched 29 Si sample, the difference goes down to one 
order of magnitude (3rd row); while for an isotopically 
enriched 28 Si sample (4th row), the difference goes up to 
three orders of magnitude. The much smaller possible 
random Overhauser field in Si QDs leads to much longer 
inhomogeneous broadening time T 2 *, ranging from 110 ns 
for a pure 29 Si dot to 11 ps for a 99.99% purified 28 Si 
dot. The long T 2 time indicates that in Si coherent elec- 
tron spin manipulation is possible without spin echoes. 
Conversely, if an inhomogeneous field is required to ma- 
nipulate the electron spin states^ ' 12 i 40 in Si such a field 
probably needs to be applied. 

We also estimate the anisotropic HFI of a nuclear spin 
at Rj in the QD by including two contributions. One is 
the near-field contribution, dominated by the anisotropy 
of the spin density for r ss Rj . Taking a small volume 
fi ~ Qpc around Rj, containing Nq ~ 2 Si atoms, the 
near-field contribution is fl \F(Rj)\ 2 Ojv n , and the scaling 
bpf = (1.1/JV) mT can be used. The second is the far-field 
contribution due to nuclear spins randomly located inside 
the QD. Our numerical estimates indicate that, for a QD 
with a radius of 20 nm, the near-field contribution is in 
the order of 10~ 4 Gauss, about 3% of the contact HFI, 
as expected. The far-field contribution is even smaller, 
at only ~ 10 -6 Gauss per nuclear spin. In short, in a 
Si QD where a single electron has a smooth probability 
distribution, the anisotropic HFI is negligibly small. 



V. CONCLUSIONS 

In summary, we have performed a state-of-the-art all- 
electron calculation of hyperfine interaction in the Si con- 
duction band. Our study introduces and validates an 
ab-initio approach to an open theoretical question of key 
relevance in spin behavior in a variety of materials, in 
particular semiconductors, paving the way for such cal- 



culations in other crystalline systems. Our calculated 
HFI strengths are consistent with existing experimental 
observations. The theoretical estimate for electron spin 
T 2 * dephasing time in a natural Si QD is one to two or- 
ders of magnitude longer than that in a GaAs QD. The 
corresponding electron spin decoherence time should be 
at least two orders of magnitude longer, formally and 
quantitatively verifying the advantages of Si as a host 
material for spin qubits. 



Acknowledgments 



This work was partially supported in Brazil by CNPq, 
CAPES, FAPESP, and FAPERJ. BK, RBC, and HMP 
performed this work as part of the INCT Program. XH 
acknowledges support by NSA/LPS through ARO with 
grant number W911NF-09- 1-0393. SDS acknowledges 
support by NSA/LPS. 



1 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

2 S. Das Sarma, R. de Sousa, X. Hu, and B. Koiller, Solid 
State Commun. 133, 737 (2004). 

3 B.E. Kane, Nature 393, 133 (1998). 

4 A.M. Tyryshkin, J.J.L. Morton, S.C. Benjamin, A. Arda- 
van, GAD. Briggs, J.W. Ager, and SA. Lyon, J. Phys.: 
Condens. Matter 18, S783 (2006). 

5 E. Abe, A.M. Tyryshkin, S. Tojo, J. J. L. Morton, W. M. 
Witzel, A. Fujimoto, J. W. Ager, E. E. Haller, J. Isoya, S. 
A. Lyon, M. L. W. Thewalt, and K. M. Itoh, Phys. Rev. 
B 82, 121201 (2010). 

6 M. Fuechsle, S. Mahapatra, FA. Zwanenburg, M. Friesen, 
MA. Eriksson, M.Y. Simmons, Nature Nanotechnology 5, 
502 (2010). 

7 A. Morello, J.J. Pla, FA. Zwanenburg, K.W. Chan, K.Y. 
Tan, H. Huebl, M. Mottonen, CD. Nugroho, C.Y. Yang, 
J. A. van Donkelaar, A.D.C. Alves, D.N. Jamieson, C.C. 
Escott, L.C.L. Hollenberg, R.G. Clark RG, and A.S. Dzu- 
rak, Nature 467, 687, (2010). 

8 B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 
027903 (2001). 

9 M. Mitic, KD. Petersson, M.C. Cassidy, R.P. Starrett, 
E. Gauja, A.J. Ferguson, C. Yang, D.N. Jamieson, R.G. 
Clark, and A.S. Dzurak, Nanotechnology 19, 265201 
(2008); 

10 M.Y. Simmons, F.J. Ruess, K.E.J. Goh, T. Hallam, SR. 
Schofield, L. Oberbeck, N.J. Curson, A.R. Hamilton, M.J. 
Butcher, R.G. Clark, and T.C.G. Reusch, Molecular Sim- 
ulation 31, 505 (2005). 

11 A. Fuhrer, M. Fuchsle, T.C.G. Reusch, B. Weber, and 
M.Y. Simmons, Nano. Lett. 9, 707 (2009). 

12 R. Hanson, L.P. Kouwenhoven, JR. Petta, S. Tarucha, and 
L.M.K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2008). 

13 M. Friesen, P. Rugheimer, D.E. Savage, M.G. Lagally, 
D.W. van der Weide, R. Joynt, and MA. Eriksson, Phys. 
Rev. B 67, 121301 (2003). 

14 H.W. Liu, T. Fujisawa, Y. Ono, H. Inokawa, A. Fujiwara, 
K. Takashina, and Y. Hirayama, Phys. Rev. B 77, 073310 
(2008). 

15 E.P. Nordberg, GA. Ten Eyck, H.L. Stalford, R.P. Muller, 
R.W. Young, K. Eng, LA. Tracy, KD. Childs, JR. Wendt, 
R.K. Grubbs , J. Stevens, M.P . Lilly, MA. Eriksson, and 
M.S. Carroll. larXiv:0906.3748l To appear in Phys. Rev. B. 

16 C.B. Simmons, M. Thalakulam, N. Shaji, L.J. Klein, H. 
Qin, R.H. Blick, D.E. Savage, M.G. Lagally, S.N. Copper- 
smith, and MA. Eriksson, Appl. Phys. Lett. 91, 213103 
(2008). 

17 W. M. Witzel, R. de Sousa, S. Das Sarma, Phys. Rev. B 
72, 161306(R) (2005); W.M. Witzel and S. Das Sarma, 



ibid. 74, 035322 (2006). 

18 W.M. Witzel, X. Hu, and S. Das Sarma, Phys. Rev. B 76, 
035212 (2007). 

19 J.L. Ivey and R. L. Mieher, Phys. Rev. B 11, 849 (1975). 

20 S.H. Park, R. Rahman, G. Klimeck, and L.C.L. Hollenberg, 
Phys. Rev. Lett. 103, 106802 (2009). 

21 W. M. Witzel, M. S. Carroll, A. Morello, L. Cywinski, and 
S. Das Sarma, Phys. Rev. Lett. 105, 187602 (2010). 

22 Ultimately, donor electron spin coherence is limited by 
electron spin dipolar interaction between donors.— In 
other words, if P donor electron spins are used as qubits 
in a Si quantum computer, it is the qubit-qubit magnetic 
dipolar coupling, or qubit cross talk, that sets the ultimate 
limit to the overall quantum coherence of the system. 

23 R.G. Shulman and B. J. Wyluda, Phys. Rev. 103, 1127 
(1956). 

24 D.K. Wilson, Phys. Rev. 134, A265 (1964). See Refs. 31 
and 32 therein. 

25 V. Dyakonov and G. Denninger, Phys. Rev. B 46, 5008(R) 
(1992). 

26 L. V. C. Assali, W. V. M. Machado, and J. F. Justo, Phys. 
Rev. B 69, 155212 (2004). 

27 R. Larico, J. F. Justo, W. V. M. Machado, and L. V. C. 
Assali, Phys. Rev B 79, 115202 (2009). 

28 R. Larico, L.V.C. Assali, W.V.M. Machado and J.F. Justo, 
Appl. Phys. Lett. 84, 720 (2004); F. Ayres, L. V. C. Assali, 
W. V. M. Machado, and J. F. Justo, ibid 88, 011918 (2006). 

29 W. Kohn, in Solid State Physics, vol. 5, ed. by F. Seitz and 
D. Turnbull (Academic Press, New York, 1957). 

30 C.G. Van de Walle and RE. Blochl, Phys. Rev. B 47, 4244 
(1993). 

31 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

32 D. J. Singh, Plane-waves, Pseudopotentials and the LAPW 
method (Kluver Academic, Norwell, 1994). 

33 P. Blaha, K. Schwarz, G.K.H. Madsen, D. Kvasnicka, and 
J. Luitz, WIEN2K, An Augmented Plane Wave + Local 
Orbitals Program for Calculating Crystal Properties (Karl- 
heinz Schwarz, Techn. Universitat Wien, Austria), 2001 
(ISBN 3-9501031-1-2). 

34 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

35 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 
(1976). 

36 Z. Wu and R.E. Cohen, Phys. Rev. B 73, 235116 (2006). 

37 C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 
(1999). 

38 B. Koiller, R.B. Capaz, X. Hu, and S. Das Sarma, Phys. 
Rev. B 70, 115207 (2004). 

39 D. Paget, G. Lampel, B. Sapoval, and V.I. Safrov, Phys. 
Rev. B 15, 5780 (1977). 



J.R. Petta, A.C. Johnson, J.M. Taylor, E.A. Laird, A. Ya- Gossard, Science 309, 2180 (2005). 

coby, M.D. Lukin, CM. Marcus, M.P. Hanson, and A.C. 



